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We study the compatibility of AMS-02 data on the cosmic-ray (CR) positron fraction with data on 
the CR electron and positron spectra provided by PAMELA and Fermi-LAT. We do that in terms of 
a novel propagation model in which sources are distributed in spiral arm patterns in agreement with 
astrophysical observations. While former interpretations assumed an unrealistically steep injection 
spectrum for astrophysical background electrons, the enhanced energy losses experienced by CR 
leptons due to the larger average source distance from Earth allow us to reproduce the data with 
harder injection spectra as expected in a shock acceleration scenario. Moreover, we show that in 
this approach, and accounting for AMS-02 results, the contribution of nearby accelerators to the 
fluxes at very high energy can be significantly reduced, thus avoiding any tension with anisotropy 
upper limits. 







Introduction: The detection of a rising positron frae- 
tion (PF) above 10 GeV, as measured by PAMELA Q, 
Fermi-LAT ji] and AMS-02 0], is one of the most strik- 
ing recent results in astroparticle physics. In the stan- 
dard scenario, electrons are mainly accelerated in cosmic 
ray (CR) sources, while positrons are only produced by 
collisions of CR protons and Helium in the interstellar 
medium. The fraction of those secondary e+ to the total 
e~ -I- e+ (CRE) flux reaching the solar system should de- 
crease with energy because of the decreasing escape time 
of the primary nuclei from the Galaxy as inferred by the 
observed secondary/primary nuclei ratios. 

Fermi-LAT observations of a rising e+ absolute spec- 
trum [J| , prove the presence of an excess in the e+ chan- 
nel and, together with the CRE spectra [5|, |6| and the 
e~ spectrum [2| measured by Fermi-LAT and PAMELA 
respectively, exclude that the rise of the PF is due to a 
steep e~ spectrum. 

The presence of an extra primary CR positron com- 
ponent at very high energy either of astrophysical or 
of exotic origin, seems therefore unavoidable. Follow- 
ing this approach, several analyses succeeded in consis- 
tently reproducing the PF measured by PAMELA, the 
Fermi-LAT CRE spectrum and other independent CR 
data sets (see Q for a comprehensive review). Most of 
them were based on the standard assumption, common 
to both semi-analytical and numerical models, that CR 
sources are distributed smoothly in the Galaxy. 

An important problem arising in this approach is that 
the required injection spectrum of the primary e~ com- 
ponent is very steep: If a power-law spectrum is as- 
sumed according to standard acceleration theory, the 
spectral index lies in the range 2.6 -^ 2.7 depending on 



the details of the propagation model. This is signifi- 
cantly steeper than that inferred from radio observations 
of SNRs, (7) = 2.0 ±0.3 [I]. Moreover, these values for 
the slope are quite different from the values 2.2 -f- 2.4 re- 
quired to reproduce the CR nuclei spectra. Therefore, it 
is very difhcult to reconcile this scenario with shock accel- 
eration theory which generally predicts the same spectral 
index, close to 2 -^ 2.3, for electrons and nuclei [0]. 

However, the Galaxy has a spiral arm structure where 
also astrophysical CR sources are more likely to be found, 
while the solar system lies in an under-dense region be- 
tween two spiral arms. Therefore, high energy CREs 
experience more energy losses than estimated within the 
standard scenario, due to the increased average distance 
they have to propagate through. This naturally softens 
their observed spectrum, offerin g th us a realistic alterna- 
tive to steep injection (see e.g. [10[, although the model 
discussed there is in strong tension with the Fermi-LAT 
e+ spectrum). 

In this Letter we use a newly developed 3-D propa- 
gation code (Dragon. v3 [ll[) to account for the spiral 
arm distribution of CR astrophysical sources and prove 
the effectiveness of this mechanism. We will show, for 
the first time, that we can reproduce the observed spec- 
tra of CR species with a primary e~ injection index close 
to the one used for nuclei and in rough agreement with 
radio observations of SNRs. We will assume that the 
extra-component sources have the same spatial distribu- 
tion as the primary ones. This may be compatible with 
enhanced secondary production in aged SNRs, as sug- 
gested in [I2I, |l3[, or in pulsars 14]. 

Furthermore, in our model a prominent contribution 
from local sources is not required. We avoid then any ten- 



sion with two relevant features of the CRE data: a) the 
smoothness of the observed CRE spectra, confirmed also 
with high accuracy by AMS-02, in contrast to the bumpi- 
ness expected if several local sources contribute to the 
high-energy fluxes; b) the stringent upper limits placed 
by Fermi-LAT and AMS-02 (in this case for the PF only) 
on the CRE dipole anisotropy, which are already chal- 
lenging, albeit not excluding yet, the local source scenario 
ISl ll6| . Remarkably, AMS-02 results suggest the extra- 
component to be softer than inferred from PAMELA ob- 
servations [17H19J- This makes the relative contribution 
of local sources, if any, less relevant than previous anal- 
yses suggested. We note that, although several SNRs 
and pulsars are observed in the nearby region {d < few 
hundred pc), none, or only a few, of these sources, may 
significantly contribute to the observed e flux, either 
because they are not powerful enough or because propa- 
gation may take place along streams which have a small 
probability to intersect the solar system J20J . 

Concerning the low energy part of the CRE spectra, 
where heliospheric propagation effects are relevant, we 
use the recently developed HelioProp code |2l| which 
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solves the CR propagation equation in the solar system 
accounting for charge-dependent drifts (see also |22|). 

In this way, we describe consistently the CR fluxes at 
Earth exploiting 3D propagation models on both galactic 
and solar system scales. 

Propagation setup: The CR transport in the Galaxy 
is described by the well-known diffusion-convection- 
reacceleration-energy-loss equation [23|, [24 1. For each 



CR particle, we solve the set of coupled transport equa- 
tions with the numerical code DRAGON [III llfl l25[ , suit- 



ably extended to describe generic 3D spatial geometries. 
The code can deal with arbitrary CR source and gas 
distributions, as well as generic magnetic field models 
and fully anisotropic, position dependent diffusion. For 
our purposes, we can take the diffusion coefficient as a 
scalar with the following dependences on the rigidity p: 
Dip) = /3-"-4 Do (p/po) , with po = 3 GV 0. 

In this work we consider a single benchmark propa- 
gation setup characterized by the following parameters: 
Do = 3 X 10^* cmVs, S = 0.6, half-halo height of 4 kpc 
and no-reacceleration. For all nuclear species we use an 
unbroken power-law source spectrum with the same spec- 
tral index 7o.p ~ 2.28. For the e~ background source 
spectrum we assume a broken power-law as required to 
consistently reproduce the spectrum of the diffuse ra- 
dio emission of the Galaxy and high energy CR data 
[26|,l27|: below 4 GeV we adopt a spectral index 1.2 while 
above that energy we tune it against PAMELA data 
(see below). This model provides a very good combined 
fit of B/C and proton spectra respectively measured by 
HEAO-3 ^ and PAMELA and also matches ^He 
and other nuclear species absolute and relative spectra 
provided the solar modulation treatment described in |2l| 
is consistently used. We remark that different choices of 
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FIG. 1. Top view of the propagated distribution, normalized 
to its maximum, on the galactic plane of e~ at 100 GeV for 
sources distributed in the spiral arms (top panel) or smoothly 
(bottom panel). In the first case the contour of the assumed 
source distribution is superimposed (black lines). 



the relevant parameters would not significantly affect our 
results, provided that all the above datascts arc consis- 
tently reproduced. We verified that introducing a spec- 
tral break a,t E ^ 200 GeV in order to reproduce the 
spectral hardening of the proton and Helium spectra ob- 
served in PAMELA and CREAM J3C| measurements docs 
not have a significant effect on the positron and electron 
spectra in the energy range we consider in this work. 

With the new 3D code, we can release the azimuthal 
symmetry hypothesis (rotational invariance about the 
Galactic disk axis) and introduce a realistic spiral arm 
pattern for the source distribution. Here we adopt the 
model used in [31| , which reproduces the observed spiral 
structure of the Milky Way. We show in Fig. [T]a face-on 
view of the density of propagated CR primary e~ on the 
galactic plane at 100 GeV, together with the distribution 
of their sources, for the two cases with and without the 
spiral arm distribution. The effect of the source distri- 
bution is striking: In the spiral arm case the electrons 
are more closely attached to their parent arm, while in 
the other case they are more uniformly spread. We veri- 



fied explicitly that the spiral arm structure does not af- 
fect significantly the spectrum of protons and nuclei, nor 
that of secondary (low energy) electrons and positrons, 
as expected because of the longer mean- free-path of these 
particles with respect to high energy leptons. 

We treat solar modulation with the recently developed 
HelioProp numerical code which solves the CR trans- 
port equation in the heliosphere accounting for charge- 



dependent drifts [21| . Before they are detected at Earth, 
CRE lose energy due to the solar wind while diffusing in 
the solar system [3^ . Because of drifts in the large scale 
gradients of the solar magnetic field (SMF), the modu- 
lation effect depends on the particle charge including its 
sign |33|. Therefore, it depends on the polarity of the 
SMF, which changes periodically every ~11 years |34| . 

The SMF has also opposite polarities in the north- 
ern and southern hemispheres: at the interface between 
opposite polarity regions, a heliospheric current sheet 
(HCS) is formed (see e.g. [H). The HCS swings then 
in a region whose angular extension is described phe- 
nomenologically by the tilt angle a. The magnitude of 
a depends on solar activity. Since particles crossing the 
HCS suffer from additional drifts because of the differ- 
ent orientation of the magnetic field lines, the intensity 
of the modulation depends on the extension of the HCS. 
The geometry of drifts is such that particles having the 
same charge sign as the polarity of the magnetic field ar- 
rive to Earth preferentially from the poles, and particles 
having opposite charge sign mainly drift inwards along 
the HCS [36|. This results into different energy losses of 
particles with opposite charge signs. Besides the magni- 
tude of the tilt angle, another very important parameter 

we assume 
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of the model is related to diffusion. As in 
that diffusion occurs in the Bohm regime, and that its in- 
tensity is described as D{p) = X{p)v/i, with A being the 
momentum-dependent mean- free-path. 

Results: We first assume the extra-component sources 
to be located in the spiral arms and have the same dis- 
tribution as standard SNRs, and disregard, at least for 
the moment, the possible role of local sources. For the 
extra-component source spectrum we assume the form 
JEc{e^) oc ^-Tf.Ec exp{-E/Ecut) and tune the involved 
parameters against the data. We normalize the primary 
electron component to the PAMELA e~ fiux at 33 GeV. 
We find that in these conditions the AMS-02 PF and the 
PAMELA and Fermi-LAT e~ spectra can consistently be 
reproduced if 70, EC — 1-75 and Ecnt — 10 TeV (see Fig.[2] 
and [3]). Remarkably, passing from a smooth source dis- 
tribution to a more realistic spiral arm pattern, a harder 
e~ source spectral index is required: 7o,bkg — 2.38, to 
be contrasted with 7o,bkg — 2.65 used, e.g., in [l6|. As 
we already pointed out, this is a consequence of the so- 
lar system being placed between two main arms (Perseus 
and Sagittarius-Carina), hence in a source under-dense 
region. This turns into a larger average distance, hence 
stronger losses, between the bulk of sources in the arms 
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FIG. 2. The e~ -f e"*" (blue), e~ (purple) and e"*" (red) prop- 
agated spectra computed in the hypothesis that the back- 
ground and extra-component sources are distributed in the 
spiral arm as shown Fig. [T] and for a SNR-like contribution 
to the extra-component. Dashed lines are for the case of a 
pulsar-like contribution. PAMELA e'^ data have been de- 
rived (without error propagation) starting from the PF and 
e~ spectrum released by the same collaboration. 
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FIG. 3. The PF computed assuming the extra-component 
sources are located only in spiral arms. The blue and green 
curves correspond to the AMS-02 and AMS-01 data taking 
periods respectively. The solid curves are for a SNR-like con- 
tribution at high energy. The dashed curve is for a pulsar-like 
contribution. The dotted line is the interstellar PF. 



and the observer. The e"*" spectrum measured by AMS-01 
and that computed on the basis of PAMELA PF and e^ 
spectrum (preliminary PAMELA e+ results agree with 
this estimate) are nicely matched by our model. We also 
notice that the spectral steepening found by H.E.S.S. [37| 
is naturally reproduced with a very high energy cutoff as 



that expected in the scenario envisaged in [12|, 113 . 



where e are produced as secondaries of CR nuclei and 
re-accelerated in SNRs. Pulsars would hardly provide 
those high energies. The pulsar scenario however cannot 



be excluded on the basis of this data. Indeed, an alter- 
native model also compatible with the pulsar scenario, 
with a lower cutoff energy E^ut — 1 TeV and a slightly 
harder spectral index 70, ec — 1-5 can reproduce the data 
comparably well, as also shown in Fig. [5] and [31 In this 
latter case, the PF is expected to flatten out and decrease 
above ^ 300 GeV. Remarkably, observations of the PF 
at these high energies could distinguish the two models. 

However, the CRE spectrum measured by Fermi-LAT 
is not satisfactorily reproduced in these conditions. This 
discrepancy can point to the presence of unknown sys- 
tematics in Fermi-LAT data. However, it could also be a 
signal of the emergence of a nearby e~ source at high en- 
ergy. We show in Fig. |4] and [5] the CRE and e~ spectra 
and the PF computed for the same model used above, 
but normalized this time to the CRE flux of Fermi-LAT 
at 33 GeV (in contrast to the secondary CRE spectra, 
whose normalization is fixed by the one of the proton 
spectrum, the normalization of primary e~ is a free pa- 
rameter of the model). 

While the PF is correctly reproduced, the CRE spec- 
trum above -^200 GeV is underproduced. This leaves 
room for a possible contribution of e~ by a local source 
at high energies. This is shown as the triple-dot-dashed 
curve in Fig. 21 The injection spectrum of this com- 
ponent is assumed to be a power-law with exponential 
cutoff, with 7o,L-EC = 2.1 and i^cut,L-EC = 1 TcV. The 
required energy output is 3.6 x 10^'^ erg if the source is 
located at a distance d ~ 290 pc from the solar system, 
compatible with the position of the Vela SNR. The dipole 
anisotropy due to this single source is diluted in the sea of 
the galactic extra-component and is ^ 0.8% at 100 GeV, 
reaching a maximum of ~ 2 % at 500 GeV, which is sig- 
nificantly smaller than present experimental upper limits 
(see also [391 for a discussion of the role of anisotropy 
in view of the CTA [40[ in the local source scenario). 
At low energy, we compute our spectra using the solar 
propagation parameters appropriate for the data-taking 
period of each dataset. In particular, we use for AMS-01 
a ~ 10° and positive polarity, while for PAMELA and 
AMS-02 we use a = 10° and a = 60° respectively, and 
negative polarity. The only free parameter is then the 
normalization of the mean-free-path A, that we tune on 
the data. Remarkably, we achieve also a very good fit of 
all PF datasets in both our models (see Fig. [3] and [S|) by 
tuning A = 0.2 and A = 0.4 at 1 GeV for AMS-02 and 
both PAMELA and AMS-01, respectively. These values 
are in rough agreement with findings from the analysis 
of the time dependence of proton spectra [41| . 

Conclusions: We have computed the CR electron and 
positron spectra at Earth within a 3D numerical propa- 
gation model. For the first time, we have computed the 
consequences of the CR sources being mainly distributed 
in spiral patterns in the Galaxy, with the solar system 
lying in an inter-arm region. The enhanced energy losses 
implied by the greater average distance that e^ have to 
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FIG. 4. The e~ -f e+ (solid blue), and e~ (solid purple) prop- 
agated spectra computed assuming that the extra-component 
sources are located only in spiral arms and that one nearby 
electron accelerator is present (triple-dot-dashed curve). The 
dotted curve represents the background e^ component. 
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FIG. 5. The PF under the same hypothesis as Fig. H) The 
line notation is the same as Fig. [S] 



cross to arrive at Earth, allow us to fit the observed CRE 
spectra and positron fraction, including the new data re- 
cently provided by AMS-02, by adopting a CRE injection 
index very close to the one needed for CR protons and 
other nuclear species. The resulting picture is, for the 
first time, compatible with expectations from shock ac- 
celeration scenarios. 

At high energies, compatibility with PAMELA and 
AMS-02 observations is achieved with a charge symmet- 
ric extra-component whose sources are distributed in the 
spiral arms. SNRs as in envisaged in [1^ [isl, |38[ or pul- 
sars [ij] , may be viable source candidates. Measurements 
of the PF at '^ 500 GeV could distinguish the two possi- 
bilities. This scenario naturally avoids tension with the 
smoothness of the observed spectra and with constraints 
on dipole anisotropics. However, Fermi-LAT data are 



not well reproduced in this model. Barring unaccounted 
systematics, compatibility with Fermi-LAT data can be 
obtained by adding a contribution of only electrons as it 
may be expected from a local SNR (e.g. Vela). 

At low energies, wc used charge-dependent solar mod- 
ulation to successfully reproduce data taken in different 
periods of solar activity. This further strengthens the 
need for an accurate description of solar propagation to 
interpret data below 10 GeV. 
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